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Abstract 

We investigate the interaction of strongly correlated electrons with phonons in the frame of 
the Hubbard-Holstein model. The electron-phonon interaction is considered to be strong and is an 
important parameter of the model besides the Coulomb repulsion of electrons and band filling. This 
interaction with the nondispersive optical phonons has been transformed to the problem of mobile 
polarons by using the canonical transformation of Lang and Firsov. We discuss in particular the 
case for which the on-site Coulomb repulsion is exactly cancelled by the phonon-mediated attractive 
interaction and suggest that polarons exchanging phonon clouds can lead to polaron pairing and 
superconductivity. It is then the frequency of the collective mode of phonon clouds being larger 
than the bare frequency, which determines the superconducting transition temperature. 
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I. INTRODUCTION 



Since the discovery of high-temperature superconductivity by Bednorz and Miiller,§ the 
Hubbard model and related models such as RVB and t-J have widely been used to discuss the 
physical properties of the normal and superconducting statci'l'Si'l However, a unanimous 
explanation of the origin of the condensate in high-temperature superconductors has not 
emerged so far. One of the unsolved questions is in how far can phonons be involved in 
the formation of the superconducting state. In experimental and theoretical works mostly 
the change of phonon frequencies and phonon life times associated with the superconducting 
transition were discussed. For example, the decrease of frequencies of Raman-active phonons 
at the transitionB and observation of the isotope effect for the case of not optimally doped 
superconductors!, as well as the observation of phonon-induced structure in the tunnel 
characteristics,! speak in favor of strong electron-phonon coupling in the cuprates. 

The aim of the present paper is to gain further insight into the mutual influence of 
strong on-site Coulomb repulsion and strong electron-phonon interaction by using the single- 
band Hubbard-Holstein model and a recently developed diagrammatic approach. 
For simplicity we consider coupling to dispersionless phonons only, although this might 
not be the most interesting case with respect to superconductivity. However, previ- 
ous investigations0@0 have shown that the Hubbard-Holstein modelSlii constitutes a 
formidable problem of its own. Other authors have also intensively studied this model 
Hamiltonian.SiiH 

Because the interactions between electrons and electrons and phonons are strong, we 
include the Coulomb repulsion in the zero-order Hamiltonian and apply the canonical trans- 
formation of Lang and Firsovil in order to eliminate the linear electron-phonon interac- 
tion. In the strong electron-phonon coupling limit the resulting Hamiltonian of hopping 
polarons (i.e., hopping electrons surrounded by clouds of phonons) can lead to an attractive 
interaction among electrons being mediated by the phonons. In this limit the chemical po- 
tential, on-site Coulomb energy as well as the frequency of the collective mode of phonon 
clouds (which is much larger than the bare frequency of the Einstein oscillators) are strongly 
renormalizedS'^ii affecting the dynamical properties of the polarons and the character of 
the superconducting transition. This will be discussed by assuming that renormalized on- 
site Coulomb repulsion and attractive electron-electron interaction completely cancel each 
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other. We suggest that the resulting superconducting state with polaronic Cooper pairs is 
mediated by the exchange of phonon clouds during the hopping processes of the electrons. 



II. THEORETICAL APPROACH 

A. The Lang-Firsov transformation of the Hubbard-Holstein model 

The initial Hamiltonian of correlated electrons coupled to optical phonons with bare 
frequency uq is given by 

H = He + Hpf^ + He-ph, (1) 

K = ^ {t{j -i) -eo 5ij} a]^ai„ + n^rii^, (2) 

ija i 

'^Ph = Yl {^i^' + i) ' '^^e-ph = g^niqi, (3) 

i i 

Ui = ^rii^, ni„ = al^ai„, qi = ^ (bi + . (4) 

Here a-^ (ctjo-) and 6/ are creation (annihilation) operators of electrons and phonons, 
respectively; i refers to the lattice site and a to the spin; is the phonon coordinate and 
g the electron-phonon interaction constant; cq = cq — /i with local energy eo and chemical 
potential /i; U the on-site Coulomb repulsion; t{j — i) is the two-center transfer integral. 
The Fourier representation of t{j — i) is connected to the tight-binding dispersion £:(k) of 
the bare electrons, 

- ^) = ^ E ^(k) exp{-zk(R, - Ri)}, 

k 

with band width W . The energy scale of this model is fixed by the parameters W , U, g and 
hujo- As an additional parameter we have the band filling. 
After applying the Lang-Firsov transformation^ 

with 
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where pi is the phonon momentum and g the dimensionless interaction constant, the polaron 
Hamiltonian is obtained as 



nj 'lyO I 'lyO I nj 

rip — rLp + rLpf^ + rtint, 

= ^ Hip, H^p = e'^rii^ + U n^Uii , 



(7) 



(8) 



a 



Hint = X] t{j - i) 




(9) 





(10) 



e = 60 — 11 = 11 + ahujQ, U = U — 2ahu!o, a = ^g"^. 



(11) 



In order to derive the polaron Hamiltonian it was necessary to include the shift of the 
phonon coordinate q- of the form e^qie~^ = Q^i which is responsible for the elimination 
of the linear electron-phonon interaction. The polaron Hamiltonian is by nature a polaron- 
phonon operator, i.e., the creation operator cj^^ and destruction operator Qo- in Tip should 
be interpreted as creation and destruction operators of polarons (electrons dressed with 
the displacements of the ions) which couple dynamically to the momentum of the optical 
phonon. In zero-order approximation (omitting T-Lint) polarons and phonons are localized 
with strongly renormalized chemical potential Jx and on-site Coulomb interaction U . The 
operator Hint describes the tunneling of polarons between the lattice sites, i.e., tunneling of 
electrons surrounded by clouds of phonons. 

B. Expansion about the atomic limit 

The problem is now to deal properly with the impact of electronic correlations on the 
polaron problem. This can be done best by using Green's functions provided one finds 
a key to deal with the spin and charge degrees of freedom. For the general case when 
U is different from zero, the Coulomb interaction has to be included in the zero-order 
Hamiltonian. As a consequence conventional perturbation theory of quantum statistical 
mechanics is not an adequate tool because it relies on the expansion of the partition function 
about the noninteracting state (use of traditional Wick's theorem and conventional Feynman 
diagrams). Similar is the situation for composite particles like polarons, Qo- = Ojo- exp(i^pj), 
involving operators for the electron and phonon subsystem. 
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It was Hubbardi who proposed a graphical expansion for correlated electrons about the 
atomic limit in powers of the hopping integrals. This diagrammatic approach was reformu- 
lated in a systematic way for the single-band Hubbard model by Slobodyan and Stasyuk,!! 
and, independently, by Zaitsevil and further developed by lzyumov.0 In these approaches 
the complicated algebraical structure of projection or Hubbard operators was used. There- 
fore, it appeared to be more appropriate to develop a diagrammatic technique which uses 
more simple creation and annihilation operators for electrons at all intermediate stages of 
the theory (for details see Ref. 10,11). In the latter approach the averages of chronologi- 



cal products of interactions are reduced to n-particle Matsubara Green's functions of the 
atomic system. These functions can be f actor ized into independent local averages by using a 
generalization of Wick's theorem (GWT), which takes into account strong local correlations 



(details are given in Ref. 10,11, |25|) . The application of GWT yields new irreducible on-site 
many-particle Green's functions or Kubo cumulants. These new functions contain all local 
spin and charge fluctuations. The analogical linked-cluster expansion for the Hubbard model 
around the atomic limit was recently reformulated by Metzner.0 

C. Averages of phonon operators 

We define the temperature Green's function for the polarons in (7) in the interaction 
representation by 

e;(x,a,r|x',a',r') = -{T c^,{t)c^,,,{t')U{I3)%, (12) 

with 

with = Tip + Tipf^ and evolution operator given by 

U{(3) = Texp jyrHWr)^ , (13) 

and where x, x' are the site indices and r, r' stand for the imaginary time with < r < /5; T 
is the time ordering operator and f3 the inverse temperature. The statistical average {■ ■ ■)q 
is evaluated with respect to the zero-order density matrix of the grand-canonical ensemble 
of the localized polarons and phonons 

11 rr„„-/3WL rr^^.-W?.. 



Tr e-^^" J- .J- Tr 6''^'^'^^ Tr e-'^'^'^p^ 
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The upper index 'c' in (13) means that only connected diagrams must be taken into ac- 
count. The density matrix (14) is factorized with respect to the lattice sites. The phonon 
part is easily diagonalized by using the free phonon operators, bi and b], while the on-site 
polaron Hamiltonian contains the polaron-polaron interaction which is proportional to the 
renormalized parameter U, which only can be diagonalized by using Hubbard operators.^ 
At this stage no special assumption is made about the quantity U and its sign and we will 
set up the equations of motion for the dynamical quantities for this general case. However, 
a detailed inspection of the equations will only be undertaken for the special case U = 0. 

Wick's theorem of weak-coupling quantum field theory can be used when evaluating 
statistical averages of phonon operators like, for example, the propagator for the phonon 
cloud, 

<l>(ri|r2) = <l)(ri - T2) = {T exp{ig [pin) - p{t2)]})o 

= exp{-lg^{T[pin)-p{r2)]\) 

= exp{-a{P) + ai\n-T2\)), (15) 



'^{n,T2\Ts,T4) = { T exp{ig [pin) + p{t2) - ^(rg) -p{n)]})o = 

= exp {-lg^{T[p{n)+p{T2)-p{T,)-pin)f)o) 

= exp{(7(|ri - Tgl) + a{\Ti - r^l) + a{\T2 - Tg]) 

+ cri\r2 - n\) - o-(|ri - raj) - cr{\r3 - n\) 

- 2a{P)}, (16) 



where 



f^dn-TsI) = g {Tp{n)p{T2))o 



cosh ( hu!o { ^ — Wl — T2\ 



. f phujo ^ 

smh 
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We now have to discuss the problem of how to calculate chronological averages of com- 
binations of polaron operators. Here we will make use of the above mentioned new diagram 
technique and the GWT. This approach has many-particle on-site irreducible Green's 
functions as main element of the diagrams. 



III. POLARON AND PHONON GREEN'S FUNCTIONS 



In zero-order approximation the one-polaron Green's function is of the form 

Gpix\x') = -(rCxa(T)CxV'(7-'))o 

= -{Ta^^{T)a^,^,{r'))o^r\T') 

= e;(°'(a;|x')$(r|r'), (18) 

with X = (x, (T, r). The simplest new element of the diagram technique is the two-particle 
irreducible Green's function or Kubo cumulant, which is equal to 

02 ''{XI,X2\X3,X4) 

X ^2°^*''(o-i,ri;o-2,r2|o-3,r3;o-4,r4), (19) 

where 

= (Tc^^(ri) c^2(r2)c^3(r3) c^^{n))o 

+ {Tc^^{ti)c^^{ts))o{Tc^^{t2)c^^{t4))o. (20) 

The first term on the right-hand side of Eq. (21) is of the form 

(Tc^,(ri) c^2(r2) c^-^in) Ca^in))o 
= (Ta^,(ri)a^2(r2)a^3(r3)a^4(r4))o$(rir2 1x3X4) (21) 

When the number of polaron operators increases more complicated irreducible Green's 
Functions like Qn^^^{xi . . . Xn\x'i . . . x'^) with n > 3 and all possible terms of their products 
will appear. The sum of all strongly connected diagrams (i.e., those which cannot be divided 
into two parts by cutting a single hopping line) containing all kinds of irreducible Green's 
functions in the perturbation expansion of the evolution operator, defines the special function 
Z{x\x') (for details see Ref. p!o|JTT|) . This function contains all contributions from charge and 
spin fluctuations. It allows us, together with the mass operator which is in our case the 
hopping matrix element, to formulate a Dyson- type of equation for the one-polaron Green's 
function,§§0'i0 

gix\x') = Aix\x') + A(x|l)t(l - 2)g{2\x'), (22) 
1,2 
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where 



A{x\x') = gl,^\x\x') + Z{x\x'), (23) 
t{x-x') = (5,,,/(5(r-r')t(x-x'). (24) 

Here x stands again for x, a, t and the sum is over the discrete indices including the inte- 
gration over r. Using the Fourier representation for these quantities, 



k LU„ 

^(^1^) = ^E;^E^"''™"^'^(k|.c^„), 

k LJn 
1 f+'^ 
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h/3 

A(x|za;„) = ^ / rfre^"""A.(x|r), (25) 

gives us the following form of the Dyson equation for the renormalized one-polaron Green's 
function, 

QM^u^n) = ll:^^].,. . (26) 

1 - £:(k) A(k|za;„) 

Here Un stands for the odd Matsubara frequency, ojn = (2n + 1) 7r//3. 

For further discussion of g(j{k\iuj) we need the Fourier representation of the zero-order 
one-polaron Green's function Q^"* defined in (18). In order to facilitate the investigation 
we have evaluated the propagator of the phonon cloud (16) in the strong-coupling limit 
a > l,li0'i 

$(r) = i J]e-^^"Wi'(zn„), (27) 
$(ifi„) = / rfre^^""+"(l"l) (28) 

2 J-f3 

with iln = 2mT//3. In order to find <l>(if2„) we use the Laplace approximation^ for the 
integral (28) which contains an exponential function with the parameter a. In the strong- 
coupling limit a ^ 1 we have obtained 

$(2fin) ^ ^"^^ = haujo = . (29) 
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This term is the harmonic propagator of the collective mode of phonons belonging to the 
polaron clouds. There are further terms describing anharmonic deviations. For a ^ 1 these 
terms will be omitted since they are smaller compared to the harmonic contribution. We 
obtain then for the Fourier representation of the phonon correlation function by making use 
of the Laplace approximationil and 



r2i ...Q4 

rl3 rP 



/?4 

Jo Jo 

the result 

^iQi,iQ2\i^3,i^4) ^ [Snun-M2,ni + ^n^^nJn^Mi] Hi^i)Hi^2), (32) 
which corresponds to 

Hn,T2\T3,u) ~ $(ri|r3) <^{t2\u) + $(ri|r4) ^(ralrg). (33) 

This means that we keep in the following only the free collective oscillations of the phonon 
clouds (29) which surround the polarons and use the Hartree-Fock approximation (32) and 
(33) for their two-particle correlation functions. In particular we will investigate the influence 
of absorption and emission of this collective mode (by the polarons) on the superconducting 
phase transition. 

By using the harmonic mode (29) the Fourier representation of the local polaron Green's 
function 

Oi'J (^^n) = I f_dr e^--a(°j (r) (34) 



is equal to 



+ 



+ 



iUn + Eq - - Uc 

e-PE. + iV(cUe)(e-^^« + e-^^-) 

iUn + Eq- E^+Uc 

g-^E-. ^ jV(cu,)(e"^^- + e-^^^) 



^ e-^^^+iV(cu,)(e-^^- + e-^^^) 

lUOn + E_^-E2 + UJ^ ' ' 



9 



where 



Zo = 1 + e-^^^ + e-^^- + e-^^^ , (36a) 
£;o = 0, E±^^e, E2^U + 2e, (36b) 

n(e) = (e^^ + l) VV(cUe) = (e'^"^ - l) ■ (37) 

Equation (35) shows that the on-site transition energies of the polarons are changed by the 
collective-mode energy ±ujc of the phonon clouds. The delocalization of the polarons due to 
their hopping between the lattice sites causes the broadening of the polaron energy levels. 
Equation (35) can be further simphfied for the case of small on-site interaction energy U of 
polarons. For C7 = we obtain 

r{0) (■ , I A N{u,) + 1 - n(e) Nju^) + n(e) 

{iujn — e) coth(/9a;c/2) -|- a;ctanh(/3e/2) 
{iujn - e)2 - a;2 

This function has the antisymmetric property 



(38) 



g^'jH^n\-e) = -g^'J{zu^\e) (39) 
which holds also for the renormalized polaron quantities, 

A^(-k, -io^nl - e) = -A„{k,iLJn\e), 

Gai-K -ii^nl - e) = -SaiK il^nl^)- (40) 

When taking f/ ^ we assume that the strong on-site Coulomb repulsion of polarons can be 
canceled by the attraction induced by the strong electron-phonon interaction. We consider 
this as a model case which allows to discuss in a transparent manner the polarons exchanging 
phonon clouds during hopping between the lattice sites. 



IV. TWO-PARTICLE IRREDUCIBLE CORRELATION FUNCTIONS 

In the following we discuss the influence of strong electron-phonon interaction on the 
two-particle irreducible Green's function. For U = the electronic correlation function in 
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(22) is equal to 

(T a^^ (ti) ao-2 ('^2) in) 0^4 (74))o 
= (T {ri)aa4{T4))o{T a^^ {r2)aas (t3))o 

- {T aa^{Ti)aas{T3))o{T aa2{T2)aa^{T4))o (41) 

because usual Wick's theorem works now. Making use of (33) we obtain for the two-particle 
irreducible Green's function (21) the relation, 

y2 (CTl, Ti; 0-2, r2|Cr3, Ts; (74T4j 

= S,,,aA,,as ^^SHri - n) 6?Sn^2 - rs) $(ti - T3) $(t2 - T4) 

- <^.i,.3<^..,.4 e^ff (Tl - r3) ^SHt2 - U) $(ri - r4) ^(r^ - n). (42) 

In the absence of exchange of phonon clouds by the polarons this quantity has to be zero. 
Indeed, if during the time of propagation of two polarons the electrons keep their initial 
phonon clouds, then the irreducible two-polaron Green's function (21) will vanish for the 
case U — 0. But because two electrons can be exchanged (independently of the exchange of 
phonon clouds) we obtain new contributions corresponding to two polarons with exchanged 
phonon clouds. Alternatively we may say that, in the case U = 0, Wick's theorem works 
separately for free electrons and free phonons; however, it does not apply for polarons as 
composite particles. Hence their cumulants do not vanish. 
The Fourier representation of (42), 

gi^^'^'iai, ioji, (72, ioj2\cr3, iuz, (74, iu^) 

= f ... /Vi...rfr4^f *'(ai,ri;a2,r2|a3,r3;a4,r4)e^"i^i+^'^^^-^'^^^-^^^^^ (43) 
Jo Jo 

has the form 

= (3 ^L0^+u;2m+'^i ^2°^"Vl> ^^1; 0-2, ^^^2^3, iuJ^] 0-4, iuJi) 
^a\,az '^0'2,<T4 ^iTl(T2 

((7i, iuJi] (72, iuJ2\(Tu «<^4; (^2, i^Jz)} , (44) 
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where 



-y 



2 



(3 ^ \%uj^ - ^] - e] \%uj^ - ^] - e] [r]2 + ^^2] [(i] + ^]l)2 + cu2] 

with f2i = co'2 — ciJs. The summation leads to 

tanh(/5e/2) [itui + ^^72 — 2e][2ti;^ — (zcui — e){iuj4 — e) — (zci;2 — e)(^^^3 — e)] 

X 



(45) 



[(i^i - e)2 - a;2] [(ia;2 - e)^ - ^c] [(^^3 - e)^ - ^c] [(^^4 - e)2 - uj^] 
2uc[iu!i + iuj2 — 2e]^ coth(/9co'c/2) pcu^ — (icui — e){iuj3 — e) — (?c<j2 — e)(ito'4 — e)] 

[(Zfil)^ - (2^e)2][(^^l - e)2 - ^2][(,^, _ ,)2 _ ^2] [(,^3 _ ,)2 _ ^2] [(,^^ _ ,)2 _ ^2] 



coth(/3cJc/2) 



+ 



(2Co'2 — e)(ic<j4 — e) + 3co'2 

1 



(icji - e)(zco'3 - e) + Scu^ 

- e)2 - ^2] [(^^3 _ g)2 _ ^2] 
2 



(46) 



[{iuj2 - e)2 - ^2] [(^^^ _ e)2 _ ^2 
The function Ao-io-z contains the contributions to the two-particles on-site Green's function 
from the different spin channels. The spin structure in Eq. (44) is due to the conversation 
law for the spins of the polarons. 

V. SUPERCONDUCTING PHASE TRANSITION 

In the following we check whether the polaronic system may have a superconducting 
instability in the absence of a direct attractive interaction for the polarons, i.e., for U = 0. 
In this case the attraction is only brought about dynamically by polarons exchanging phonon 
clouds. With respect to superconductivity we need in addition to the normal state Green's 
function (13) the anomalous propagators.ii For simplicity we limit the discussion to s- 
wave superconductivity as in previous investigations of superconducting instabilities in the 
Hubbard mo deli'0 and Hubbard-Holstein model in the strong-coupling limit, a ^ 1.ESI 

For a description of the superconducting state we need the three irreducible functions 
Ao-, Yo-,_(j and F-o-.o- which represent infinite sums of diagrams containing irreducible many- 
particle Green's functions. In order to obtain a close set of equations we will restrict ourselves 
to a class of rather simple contributions which, however, contain the most important charge. 



spin and pairing correlations; for details see Ref. This class of diagrams is obtained by 
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neglecting contributions for which the Fourier representation of the superconducting order 
parameters, l^,_o- and Y^a,a, depend on the polaron momentum k. In this approximation 
y^^_o. has to be obtained from 

X Q^^^^{a,iuj; —a, —iuj\a,iLJi; —a, —iui). (47) 
In the same approximation A^. has to be computed from 



- e{k) Y„_a{i^i) y-a,a{'i'^i)} ^2°^*''(cr, iuj; cr, iuJi\a, iuui] a, iuj) 

- £(-k) Y^_„{iuji)Y_^^„{iuji))Qf^''''{a, iu; -a, iu;i\ - a, iui; a, iu) (48) 

with 

^.(k, lu) = [1 - £(k) K{iu;)][l - £(-k) A-^i-iu)] + £(k) £(-k) Y,,_,{iu;) Y_,,,{iu;). (49) 

The corresponding equation for Y-rj,cr{'i^) can be obtained from the expression for Y^r^-ai'^^) ■ 
The last equations together with the equations for the one- and two-particle Green's 
functions determine completely the properties of the superconducting phase, provided it 
exists. In order to gain further insight into the physics contained in (47) and (48) we will 
linearize the equations in terms of the order parameter Y^-aiiou) which will determine the 
critical temperature Tg. The resulting equation for the order parameter is of the form 

1 £(k)5(-k)y;,_„(?M) 



PN^[l- s{k) K{iuji)][l - e{-k) A-^i-tuJi)] 
X Q2^^^{a,iLJ] —a, —iLj\a,iLJi] —a, —iuJi). (50) 
This equation must be solved together with the equation for h.cj{iuj) which may be approxi- 
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mated by setting the order parameters to zero giving, 
K{iuJ) = G^Jiiuj) -j^l^ 1 - £(k) A^(ia;;) ^ 

In order to determine we must solve (51) for and insert the result in (50). The 
irreducible functions in (50) and (51) can be written as 

Qf^^'^'ia, icu; a, iuji\a, icuf, a, icu) = ^2^2 {^'^li^ + ^i) tanh(/3e/2) 

'°*^('^'^^/^) [(XX, + uDiAA, + Sut) - 2u;l(A + A,)(xx, - u'^)] }, (52) 



ujc[{iuj -iuJiY - 4:UjI] 



Q'f''''{a,iu\-a,iui\ - a,iuji;a,iuj) = ~ ^2^2 {^c{x + xi){A + A;) tanh(/3e/2) 



+ coiHMix + x^rixx, - .1) } + coth(/3a;./2)((... + 3a;g) ^ ^^3^ 



]^^^^{a,iuj] —a, —iuulcTjiujf, —a, —iuji) = 

2e{2uc)^tanh{Pe/2)[icuicui + - u^] 

~ [^2 + + ^J2] [^2 + _ ^^)2] [^2 + + ^J2] [^2 + _ ^J2] 

2a;cCoth(/9(x'c/2)[iu;iu;; + 2u;c(e — u;c) — (e ~ ^c)^] 



+ 



+ (e _ ,^J2]^2 + _ ^^)2][(^ _ ^^)2 + + (2u;e) 



21 



where 



2a;c coth(/3a;c/2) [za; zo;; - 2LUc{e + a;c) - (e + uJcY] 
^ [u;^ + ie + u;,y][u;f + (e + u;,mu;-u;iy + {e+(2u;,yy ^ ^ 



X — iuj — e, A — {iu — ef — oj\ (55a) 

xi = iuJi - e, A; = {iuji - ef - lj"^. (55b) 
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For further discussion of (50) and (51) we introduce the following shorthand notation: 

5fo-(«<^) = ^(x = x'\iu) 

N ^ 1 - e{k) A^itujy ^ ' 

,.c/ • N ^ ± V g(k)g(-k) 

' iV^[l-5(k)A.(^^)][l-£(-k)A_,(-z^)] 

A^{iuj) - A_^{-iujy 

Furthermore we assume that £(k) = £(— k) holds with X]k^(^) — Z]k^^(^) — 0- Sums will 
be replaced by integrals, 

^5^= /"&Po(£), (59) 



2eV 1 kl<- 



Po(^) = Z^Wl-(^) x<( ' ' ^' (60) 

kl<f, 



where W is the band width and po a- model density of states of semielliptic form. Since 
we do not consider magnetic states here, the spin index in the paramagnetic phase can be 
omitted: 

Aa-{iu;) — A_o.(ia;) = A(ia;), (61a) 
(f)aiiuj) = (f)-ai'ii^) = 0(«<^), (61b) 
C(^'^) = 0-.(^'^) = 'l>'%^^)■ (61c) 

However, the spin index is essential for the superconducting order parameter Yfj^^^^iuj), 

Y^_^{iu}) = g^^^^Y{iu}), (62a) 
ga-a = Sa,^-Sa,i, (62b) 

where Y{iu!) is an even function of the frequency, 

Y(iu;) ^Y{-iu;). (63) 
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Finally we have to add the equation which determines the chemical potential: 

Here Np is the number of polarons and N the number of lattice sites. With (59) and (60) 
the functions (56) and (57) can be written as 

'/'M = 7T7—\ — = . (65) 



4 (l - VI - A^M) 



A(ia;) 



^(^^) = w w^^ = . (66) 



W X{iu;) W 



tLU 



with A(ia;) = (M^/2) A(ia;). In order to check whether the state to be determined from Eq. 
(51) is metallic or dielectric one has to analyze the renormalized density of states given by 

where X{E + iO"*") is the analytical continuation of X{iuj). 

VI. ANALYTICAL SOLUTIONS 

The expressions for A(ia;) and Y{iu!) simplify by using the shorthand notations (56) and 
(57) and the symmetry property (62): 

Y{iLu) = ^^(p'''{iuji) g'^^^''' {a,iuj; -a, -iuj\a,iuji; -a, -i(^i)Y {iuj^ (68) 



A(ia;) = g^°\iuj) - ^'^4>{m) [^^°^ (cr, ia;; cr, ia;/| cr, ia;;; cr, ia;) 

+ Q^^^^^{a,iLJ] —a,iLJi\ — a,iLJi] a,iLj)] . (69) 
In order to find a solution of Eq. (68) we insert (54), replace Y{iu;n) by 
Y{lUJn) = </>*"(^o) Xi^oon) Y{zo) with ;2 ~ 0, (70) 



2ujc[oJc - etanh(/3e/2) cosh(/3cc;c/2)] + cos\i{I3uJc){-ujI + + ujl){cosh{(5uJc) - 1)"^ 

['^^+(^c + 6)2][^^+(^c + 6)2] ■ 

(71) 
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and make use of Poisson's sum rule, 

\ y f{iun) = I dz (72) 

p^JK nj 27riJc ef^^ + V ^ ^ 

where C denotes the usual counterclockwise contour of the imaginary axis. We obtain then 

from Eq. (68) with help of the analytically continued function x{^) ior Z — Zq — an 
equation for the critical temperature Tc, 

x(O|e)0-(O|e) = l, (73) 



X(0|6) = f 2^e[a;e - 6tanh(/3e6/2) coth(/3,^e/2)] + ( 2 _ ^2^-2 _ ^^^^ 

This quantity is even in e and therefore only the absolute value of e = eo — /i determines 
ksTc = From (58) and (65) one can make a rough guess for the quantity 4>^'^{0): 

r(o)^(^)':^, 7 = ^ (1 + yr^WT^) , (75) 

where 7 has to satisfy: 7(— e) = 7(e). This quantity can be obtained self-consistently from 



Eq. (64) for the chemical potential. For simplicity we replace here [1 + ^^T^^JJiO'+iS)] by 
27. Then (64) can be written as 

^^EA(-»)e""°*=f. (76) 

Using (69) together with (52) and (53) we can express A(ia;„) as 

{iuJn - r)Ai{c) + ^\:Di{e) 



{iujn -eY -ujI 
, ujl[{iuJn-e)A2{e)+uJcB2{e)\ 



with unknown coefficients Ai and Bi. They can be found from Eq. (69) or more easily from 
the asymptotic behavior of (77) for |cj„| — > 00, 

... Aie + uJcBi Ai{ujI + £2) + 2euJcBx + UJIA2 



, Ar{e^ + ^eujl) + B^{ujl + ^e'u,)+ujl{'ieA2 + uj,B2) , 

H T-. TT h... {(Q) 
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If we compare this with the asymptotic behavior of the full one-polaron Green's function 
(Appendix A) by invoking the methods of moments together with the asymptotic behavior 
of g{iu!n) in (66), we obtain 

^i(e) = 1, (79a) 
Si(e) = --[Mi + e], (79b) 

OJc 

A2{e) = ^ [M2 + 2eMi + - - {W/4:f] , (79c) 

52(e) = ^ [-M3 - 3eM2 + Ml {luI - 3e^ + 3{W/4)^) + eul - + 2>e{WlAf] , (79d) 

where Mj is the i-th moment of the one-polaron Green's function. The results in (A. 5) 
for the moments in lowest order allow to evaluate Ai and Bi, see (A. 7). Ai = 1 describes 
the asymptotic freedom of the polarons. Bi = tanh(/3e/2) is identical with its value in the 
zero-order polaron Green's function of (38). The two new quantities A2 and B2 are small 
quantities being proportional to ujo/i^c— ^/cx 

Inserting now (77) into the left-hand part of Eq. (76) and performing the summation 
leads to 

l^A(iu;„)e-"°^=n(6) 

tanh(/3e/2) [tanh{pUc/2) - 1][1 - tanh\(3e/2)] 
2[1 - tanh^(/5a;c/2) tanh2(/3e/2)] 

B2(e) , 1 - tanh2(/?e/2) 
+ ^ tanh(/3a;e/2) ^ . ^ ^^^y" 



4 " '1 -tanh"(/3e/2)tanh'(/3cUc/2) 
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A2(e) + B2(e) A2(e) - B2(e) 



.cosh2[/?(cJc + e)/2] cosh^[/?(cJe - e)/2] J ' 
which, according to Eq. (76), is equal to (7/2) (iVp/iV). Because of the large collective 
frequency, f^ujc ^ 1, we may omit exponentially small quantities like exp(— /?a;c). Since we 
are interested in results for electron numbers which are close to half filling (e = 0), also 
|e| <^ LOc holds. Furthermore we will neglect contributions of the order 1/a. Then the 
equation for the chemical potential (74) is simply 

n(e) = 7np/2, rip = N^/N. (81) 

If furthermore we use 7 = 1 (free polarons) we obtain with (75) for 0*'^(O) the value 

(f)''=(0 = (82) 
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which allows to write the equation for the critical temperature Tc in the following form 

+ ujI - 2ea;etanh(/^;e/2) = {ujI - e^f{A/Wf. (83) 

In this approximation Tc depends only on the local parameters; but we expect that 
close to half filling this should give an impression which of the local quantities is most 
important for superconductivity in the strong-coupling limit of the Hubbard- Holstein model. 
For strictly half filling Eq. (83) can only be fulfilled when uo^ = W/4. This may perhaps be 
an unphysical large value for the renormalized quantity. It also shows that the specific limit 
tJ — is probably the critical value for the occurence of superconductivity in the frame of 
the Hubbard- Holstein model. It is clear that for C7 < superconductivity is possible; but in 
this case it would have to compete in energy with the energies of charge-ordered states. 

For the special case ojc = W/ 4 we obtain 

\2 To / / \2l 



{e/ucf [3- (eM) ] = (2|e|/u;e)tanh(/3e|e|/2). 



(84) 



Since (e/cjc < 3 holds (not discussed in detail) we may seek for solutions for the case |e| -C Wc 
leading to 



kBTc 



3 

W_ 
12 



12 



+ ... 



' 20 / e \2 



(85) 



In spite of the many approximations (which however are all reasonable) used the result 
for Tc is in so far remarkable as it shows that the critical temperature depends on the band 

width (corresponding to the largest cut-off energy of the model) and not on the effective 
mass of the ions, for small deviations from half filling Tc decreases and does not depend on 
the sign of e. 

For different values of cOc, 



OJc 



W/A - y, 



(86) 



with y 7^ there are only solutions for not half fiUing. In this case Eq. (84) may in this case 
be written as 

1 + K 



PcM — In 



K, — 



2\e\u!r ' 



< K < 1. 



(87) 
(88) 
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The condition k < 1 is equivalent to 



\e\+uJc<W/A (89) 

while «; > is related to the parameter 

< y > : - <^< e^^, (90) 

uj, > W/A, y<0: e^,„ < < e^,, (91) 



eL.,min = + (l/2)(W^/4)' ± {W/8W{W/Ay + 8u;l (92) 

For small y we can simplify (87) and (88), 

K ~ 2/(VF|e|) {£'[3 - e\4/Wf] + y[W/2 - Ae^W - e\4{Wf] 

+ y'[-3 + e\^/Wf - e\^/W)%} (93) 

with the following restrictions for e: 

y > : t/^ < < 3{W/Af - {b/Q)Wy + (29/27)?/^ (94) 
y<0: (iy/6)|y| + (25/27)y'<e'<3(iy/4)' + (5/6)iy|y| + (29/27)yl (95) 

Large values for Tc can be achieved for k ^ 1 and in the vicinity of half fiUing (e 7^ 0) : 

W5 , , W 

ksTc ^ —T-r ^ = — n < 1, y = ^c, (96) 

12((5-1)' ' ^ 4 ^ ^ 

but only for y < and hence uJc > W/i. 



VII. SUMMARY 



We have discussed the occurence of superconductivity for the Hubbard-Holstein model 
in the strong-coupling limit {g 1). Strong coupling leads to a renormalization of the one- 
polaron Green's function already in the local approximation. For ^ ^ 1 we find a collective 
mode for the phonon clouds estimated by evaluating integrals in the Laplace approxima- 
tion. Due to the absorption and emission of this mode by the polarons, the on-site energies 
of polarons are renormalized. Similarly the irreducible two-particle Green's are renormal- 
ized. Allowing for the exchange of polarons including their phonon clouds leads to a new 
irreducible Green's function which has been used to study spin-singlet pairing of polarons. 
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Analytical results for the superconducting phase have been obtained for the limiting case, 
for which the local repulsion of polarons is exactly canceled by their attractive interaction. 
The resulting equation for the critical temperature has been obtained by assuming a large 
collective-mode frequency and a nearly half-filled band case. The parameters which deter- 
mine Tc are uUc (> W/4), e (e = corresponds to half filling) and the band width W with Tg. 
In the strong-coupling limit we obtain a critical temperature which is of the order of W/12. 

Acknowledgments 

This work was supported by the Heisenberg- Landau Program and by the Program of 
the High Council for Science and Technical Development of Moldova. It is a pleasure to 
acknowledge discussions with Prof. Plakida and Prof. Crisan. V. A. M. would like to thank 
the Universities of Duisburg and Salerno for financial support. P. E. would like to thank the 
Bogoliubov Laboratory of Theoretical Physics, JINR, for the hospitality he received during 
his stay in Dubna. 



* 



APPENDIX A: METHOD OF MOMENTS 



By using the Heisenberg representation of the one-polaron Green' function 



6;,(x-x'|t-t') 



(Al) 



with 



Cxa(T") 




(A2b) 



(A2a) 



we can write the asimptotic expansion of the Fourier representation (25) for |a;, 



oo as 




(A3) 



M„ = ({4,, in[n...[n,c^^\ ]})h. 



(A4) 



n 
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Here the statistical average {■■■)h is defined with respect to the full density matrix of 
the grand canonical ensemble. In the simplest approximation we obtain for the first three 
moments of the Green's functions the relations, 

Ml = -(e + a;ctanh(/3e/2), (A5a) 

M2 = + col + {W/Af + 2ecUctanh(/?e/2) + cuqc^c coth(/?a;c/2), (A5b) 
M3 = -{e^ + 3e[u;2 ^ oj^ojQcoih{l3uJc/2) + [W/Af] 

+ u;ctanh(/3e/2) [3e^ + ?,{W/Af + o;^ + o;^ + 3a;ot^c coth(/3a;o/2) ]}. (A5c) 

The expressions for the moments can be used to determine the unknown coefficients, 
An{e) An{e), in the relation for h.^{iijj), see Eq. (78), by considering also the asymptotic 
behavior of gaii<^) in (66) for small values of Ao.(ia;), 

g^iiu;) = {2/W)X,{iu;) [ 1 + (X'/A) + 2{XyAf + ...]. (A6) 

We then insert the asymptotic form for Xaiicu) from (78). Comparing the corresponding 
equations fixes the coefficients A„(e) and A„(e), 

A(e) = 1, (A7a) 
Bi{e) = - — [Ml + e], ~ tanh(/3e/2), (A7b) 

LOc 

A2{€) = ^ [M2 + 2eMi + - a;,' - {W/Af] ~ — coth(/3u;/2), (A7c) 
B^ie) = [-M3 - 3eM2 + Mi {ujI - 3e^ + 3{W/A)^) + eul - + 2>e{W/Af] 



~ — tanh(/5e/2) 



— + 3 coth(/?a;o/2) 

UJc 



(A7d) 
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